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' Abstract. We present new yields of ^'Al from Wolf-Rayet stellar winds based on rotating stellar models which account well 

I for numerous observed properties of massive stars. We study the impacts on the yields of a change of initial mass, metallicity 

. and initial rotation velocity. We also consider the effects of a change of mass loss rates during the Wolf-Rayet phase. 

CIh' We show that for surface rotation velocities during the core H-buming phase matching to the observed ones, the quantity of 

JL , ^''Al ejected by a star of a given initial mass and metallicity is roughly doubled when the effects of rotation are taken into 

; I ■ account. The metallicity dependence of the yield is, on the other hand, very similar to that obtained from non-rotating models. 

' We estimate that at least about 20% to 50% (e.g. ~ 0.6 - 1.4 Mq) of the live ^''Al detected in the Milky- Way originates from 

, Wolf-Rayet stellar winds. We show the importance of a good knowledge of the present metallicity gradient and star formation 

IL" ' rate in our galaxy for modeling both the variation of the ^*A1 surface density with the galactocentric distance and the global 

. 5^ ' contribution of the Wolf-Rayet stellar winds to the present galactic mass of Al . 



5_( ' Key words. Stars: nucleosynthesis, abundances - rotation - Galaxy: chemical evolution 



1. Introduction 

The ^^Al radionuclide, live or extinct, is observed in at least 
three different types of media : 1) excesses of daughter ele- 
ments from ^^Al are observed in presolar Stardust grains (see 
the review by Clayton & Nittler 2004|l; 2) evidence for injec- 
tion of live ^^Al in the early Solar System nebula can be found 
in meteorites (see e.g. MacPherson et al. 1995 1; 3) ^*A1 is also 
detected in the galactic interstellar medium through its decay 
emission line at 1.809 MeV (Mahoney et al l 19821 Diehl et al. 
11995) . This diffuse emission is observed in the plane of our 
Galaxy. From its total intensity galactic ^*A1 masses have been 
derived in the range from 1 .5 Mq to 3 M© , depending much on 
the assumed underlying ^''Al density distribution, and in partic- 
ular, on the yet poorly determined scale height of ^^Al (see e.g. 
Fig. 4.7 in Knodlseder 1997). In this paper we shall focus our 
attention on the galactic interstellar ^^Al . 

Since this nuclide has a lifetime of roughly 10* yr, much 
shorter than the Galaxy lifetime, it must be continuously pro- 



duced by one or more nucleosynthetic source to be observed. 
The question which still remains to be solved is what is the na- 
ture of these sources. A coiTect answer to this question is not 
only interesting on its own, but will enable to use the emissiv- 
ity at 1.809 MeV as a penetrating tooP for studying regions of 
current nucleosynthetic activity in the MiUcy Way. 

From the observed distribution of the emission in the plane 
of our Galaxy it is generally inferred that the massive stars, 
either through their winds (Dearborn & Blake . 1985,, Vuissoz et 
al. l2004t or through their ejections at the time of the supernova 
event (Timmes et al. '1995' Woosley & Weaver '1995 ), are the 
main sources of ^*A1 (see also the review by Prantzos & Diehl 
1 19961 and references therein). Here we shall only address the 
contribution of the winds of Wolf-Rayet (WR) stars. 

Wolf-Rayet versus core collapse supernovae 

Although WR stars end their lives as supernovae, and 
also contribute to the interstellar medium enrichment of 
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Al through explosive nucleosynthesis, we distinguish the 
"wind" contribution from that of the "supernova explosion" 
considering that observations allow us to disentangle the two 
channels for ^^Al production. Indeed the ejection of -''Al by 
the winds of massive stars is not accompanied by the ejection 
of **'Fe, a radionuclide with a lifetime of 2.2 Myrs, whereas its 
ejection at the time of the supernova explosion is. Secondly, 
after a burst of star formation, the enrichment of the interstel- 
lar medium by the WR stellar winds begins prior to the en- 
richment by the supernovae explosions (Cervino et al. 2000 
Pluschke et al. I20()2> . This is due to the fact that WR stars 
originate from more massive (and thus shorter lived) progen- 
itors than the bulk of the core collapse supernovae (hereafter 
cc-SNae). Therefore observations of sufficiently young associ- 
ations of massive stars, as the Cygnus region (Knodlseder et 
al.|2002), allow us in principle to disentangle WR wind and 
cc-SNae contributions. 

Recently, Smith (2004a'( has reported a first detection 
of the 1.117 and 1.332 MeV gamma-ray lines attributed to 
the radioactive decay of ^''Fe from the inner Galaxy region, 
using the Reuven Ramaty High Energy Solar Spectroscopic 
Imager (RHESSI). After a more deeper analysis, he revised 
his first flux estimates downwards, to a marginal 2.6 sigma 
detection level (Smith, '2004b'). The observed flux in each 
of the ''"Fe line amounts to about 10 % of the -''Al flux, 
slightly smaller than the flux ratio predicted by Timmes et 
al. J1995 > from calculations of type II supernovae^ (hereafter 
SNII) nucleosynthesis by Woosley & Weaver (1995) . Taken 
at face, such an agreement would mean that all the galactic 
^^Al could originate from the SNII. However, as discussed 
in Prantzos (2004i, more recent models of the evolution and 
nucleosynthesis of cc-SNae such as those of Rauscher et al. 
( 1200 2 1 or those of Limongi & Chieffi (2003 1, predict much 
higher ''"Fe/^^Al line flux ratios. These models seem thus to 
leave more room for sources, such as the WR stellar winds, 
ejecting ^^Al and not ^Fe. According to Prantzos f2004), with 
these new stellar yields, the RHESSI results can be recovered 
provided that at least half of ^^Al originates from WR stellar 
winds. There are at least three other arguments supporting 
some additional contributions to that of the supernovae. Firstly, 
the new supernovae yields appear too low to allow an injection 
rate of 2 Mq of ^''Al per Myr (Prantzos l2004t . Secondly the 
observation of a bump of emissivity at 1.809 MeV in the 
direction of the Cygnus region, a region with no sign of recent 
supernova activity, also points towards the existence of sources 
such as the WR stars, that eject ^''Al even before the first 
supernovae appear (Knodlseder et al. 2002 1. Finally, from the 
point of view of stellar evolution models, ^^Al enrichment 
of the interstellar medium by WR stars is an inherent and 
unavoidable consequence of the most massive stars evolution 
if their observed properties are to be reproduced (Meynet et al. 
1997). An important contribution of WR stars to the galactic 
^^Al is therefore unavoidable. 



^ The mass of the progenitors of type II supernovae lies between 
llMg and 40Mo , excluding WR stars in the case of non-rotating 
models. 



The purpose of the present work is to provide yields from 
WR stellar winds based on recent grids of rotating stellar mod- 
els (Meynet & Maeder 2'Un^ EUU?! hereafter paper X and XI 
respectively). In these previous two papers we have investi- 
gated the effects of rotation on the formation of the WR stars 
at different metallicities. We showed that the inclusion of ro- 
tation makes the formation of WR stars easier, thus helping to 
reproduce the observed variations of the number ratio of WR 
to O-type stars as a function of the metallicity, of the ratio of 
WC to WN stars for metallicities below about Z = 0.020 (solar 
metallicity), as well as the observed variation with the metal- 
licity of the number fraction of type Ibc supernovae with re- 
spect to SNII (Prantzos & Boissier 2003 1. Furthermore, models 
with rotation are able to account for the small (although signif- 
icant) fraction of WR stars presenting at their surface both H- 
and He-burning products (see e.g. Crowther et al. 1995 1. None 
of these observed features can be fitted by non-rotating stellar 
models computed with recent mass loss rates prescriptions (see 
Sect. 2). 

Despite the fact that the rotating models well reproduce the 
mentioned above observed features as well as the observed sur- 
face enrichments (Heger & Langer '2000! Meynet & Maeder 
pOOO hereafter paper V) and the number ratio of blue to red 
supergiants in the Small Magellanic Cloud (Maeder & Meynet 
i2001.; hereafter paper VII), the rotating models have some dif- 
ficulties in accounting for the observed fraction of WC to WN 
stars at high metallicity. Typically, rotating models seem to un- 
derestimate this ratio at twice the solar metallicity by about a 
factor of two. In paper XI we attributed this remaining discrep- 
ancy as due to one (or both) of the following two causes: 

1) The incompleteness of the observed sample in high 
metallicity regions might bias the result. This seems already 
to be the case at solar metallicity where according to Massey & 
Johnson ( 1 998 J the observed WC/WN ratio is overestimated. 
The less luminous WC stars are more difficult to detect than 
the more luminous WN stars, especially in regions where the 
extinction may be high. In contrast, in the Magellanic Clouds, 
where the models well reproduce the observed ratio, the stars 
have nearly all the same metallicity, are at the same distance 
and the internal extinction is low. 

2) The mass loss rates during the post H-burning WN phase 
are underestimated. This would increase the WC phase dura- 
tion at the expense of the WN phase. Let us note that even if the 
mass loss rates would be underestimated in the post H-burning 
WN phases, this would not deeply modify the present results 
as far as ^''Al is concerned. 

Despite the above remaining difficulty, we think that the 
present models can sufficiently well account for the observed 
properties of massive stars to allow an acceptable estimate of 
their ^''Al yields. In Sect. 2 we briefly recall the physical ingre- 
dients of the models. The effects of rotation on ^^Al production 
by WR stars at various metallicities are discussed in Sect. 3. A 
small subset of the present models for the metallicities at Z = 
0.020 and 0.040 were already briefly discussed in Vuissoz et 
al. ( I2004> . An estimate of the global contribution of WR stars 
to the ^*A1 present in the Milky Way is given in Sect. 4, while 
Sect. 5 presents the main conclusions of this work. 
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2. Physical ingredients of the models 

The grid of models presented here was computed with the 
Geneva stellar evolution code. The physical ingredients are es- 
sentially the same as in paper X and XI, except for the inclusion 
of the Ne-Na and Mg-Al reaction chains and the omission of 
the effects of the wind anisotropics induced by rotation (see be- 
low). Let us briefly recall the main points relevant for ^^Al pro- 
duction: 

- The initial compositions are adapted for the different metal- 
licities considered here. For a given metallicity Z (in mass 
fraction), the initial helium mass fraction Y is given by the 
relation Y - Yp + AY/AZ ■ Z, where Yp is the primordial 
helium abundance and AY/AZ the slope of the helium-to- 
metal enrichment law. We use the same values as in paper 
VII i.e. Yp = 0.23 and AY/AZ = 2.5. For the metaUicities 
Z = 0.004, 0.008 and 0.040 considered in this work, we 
thus have X = 0.757, 0.744, 0.640 and Y = 0.239, 0.248, 
0.320 respectively. For the mixture of the heavy elements, 
we adopt the same mixture as the one used to compute 
the opacity tables for solar composition (Iglesias & Rogers 
[T996ll . 

- Since mass loss rates are a key ingredient of the models 
in the mass range considered here, let us recall the pre- 
scriptions used. The changes of the mass loss rates M with 
rotation are taken into account as explained in Maeder & 
Meynet ('2000a hereafter paper VI). As reference, we adopt 
the mass loss rates of Vink et al. (2000; 2001 1 who account 
for the occurrence of bi-stability limits which change the 
wind properties and mass loss rates. For the domain not 
covered by these authors we use the results by de Jager et 
al. (Q988 ). For the non-rotating models, since the empirical 
values for the mass loss rates are based on stars covering the 
whole range of rotational velocities, we must apply a reduc- 
tion factor to the empirical rates to make them correspond 
to the non rotating case. We used a value of 0.85, as in pa- 
per VII. During the WR phase we use the mass loss rates 
by Nugis & Lamers (2000). These rates, which account for 
the clumping effects in the winds, are smaller by a factor 2- 
3 compared to the ones used in our previous non-rotating 
"enhanced mass loss rate" stellar grids with predictions of 
26 Al yields (Meynet et al. [T997b . 

- During the non-WR phases of the present models, we as- 
sumed that the mass loss rates depend on the initial metal- 
licity as M(Z) = (Z/ZQ^'^MiZs) (Kudritzki & Puis ,2000, 
Vink et al. l2001l l. while during the WR phase we assumed 
no metallicity dependence. At Z = 0.040, we also com- 
puted a series of models with metallicity dependent mass 
loss rates during the WR phase. According to Crowther et 
al. (^002 ), mass loss rates during the WR phase may show 
the same metallicity dependence as the winds of O-type 
stars, i.e. scale with ~ (Z/Zo)'^^. 

- A moderate overshooting is included in the present rotat- 
ing and non-rotating models. The radii of the convective 
cores are increased with respect to their values obtained 
by the Schwarzschild criterion by a quantity equal to 0.1 
Hp, where Hp is the pressure scale height estimated at the 
Schwarzschild boundary. 



- The effect of rotation on the transport of the chemical 
species and of the angular momentum are accounted for 
as in papers VII and VIII (Meynet & Maeder 2002). All 
the models were computed up to the end of the He-burning 
phase. 

- As initial rotation, we have considered a value equal to 300 
km.s ' on the ZAMS for all the initial masses and metaUic- 
ities considered. At solar metallicity, this initial value pro- 
duces time averaged equatorial velocities on the main se- 
quence (hereafter MS) well in the observed range, i.e. be- 
tween 200 and 250 km.s ' . At low metaUicities this ini- 
tial rotational velocity corresponds also to mean values be- 
tween 200 and 250 km.s"' on the MS, while at twice the 
solar metallicity, the mean velocity is lower, between 160 
and 230 km.s ' (see Table 1 in paper XI). Presently we 
do not know the distributions of the rotational velocities 
at these non-solar metaUicities and thus we do not know if 
the adopted initial velocity corresponds to the average ob- 
served values. It may be that at lower metaUicities, the ini- 
tial velocity distribution contains a larger number of high 
initial velocities (Maeder et al. 1999 1, in which case the ef- 
fects of rotation described below would be underestimated 
at low metallicity. 

- The Ne-Na and Mg-Al nuclear reaction chains have been 
included in addition to the usual nuclear reactions for 
the H- and He-burning phases. The long-lived (fi/2 - 
7.05 X 10^ yr) ^''AF ground state is considered as a separate 
species from its short-lived (fi/2 = 6.35 s ) ^^Al'" isomeric 
state. In the following, for purpose of simplicity, the long- 
lived 26 Al* ground state will be denoted by the symbol 
26 Al . The nuclear reaction rates are taken from the NACRE 
compilation (Angulo et al. I1999I I. Arnould et al. (Il999t 
have studied in a simple parametric model the uncertain- 
ties on the abundances of some elements due to the un- 
certainties of the nuclear reaction rates. For what concerns 
26 Al they reach the interesting conclusion that the large un- 
certainties on the 26a1(/?, y)2^Si reaction rate scarcely affect 
the abundance of 26 Al , considering that "even the highest 
NACRE proton capture rates are not fast enough for lead- 
ing to a substantial destruction of the two Al isotopes by the 
time H is consumed^ (Arnould et al. 1 19991 1 . 

Let us finally recall that the wind anisotropies induced by rota- 
tion were neglected. This last choice appears justified in view 
of the results obtained in paper X. Indeed for the initial veloc- 
ities considered (wini - 300 km.s"' ), the effects of the wind 
anisotropies have been shown to be very small. Let us how- 
ever emphasize that this is not true for higher initial velocities 
(Maeder l2002b . 

3. Effects of rotation on ^^AI production by WR 
stars 

26 Al is mainly produced by proton capture on 25Mg seeds 
in the convective core during the central Hydrogen Burning 
phase (see Fig.[nand Fig.|2}. In non-rotating stars, part of the 
26 Al produced in the core, depending on the initial stellar mass, 
may appear at the surface due to the removal of the envelope by 
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Fig. 1. Evolution as a function of time of surface (^^Alj) and 
central (^^Al^) mass fraction of ^^Al for both rotating (upper 
panel) and non-rotating (lower panel) 120 Mq stellar models 
at Z = 0.020. Light dashed lines show the scaled values of the 
mass fractions 60 + ^''Al, x 2 10^ and 10 + ^''Al,. x 10^ Bold 
solid lines mimic a Kippenhahn like diagram, and show the 
evolution of the mass of the convective core and that of the total 
stellar mass as a function of time. The down-going arrows indi- 
cate the beginning of the strong wind phase, and the up-going 
arrows indicate the beginning of the surface ^''Al enrichment. 



stellar winds. From this point on, stellar winds enrich the inter- 
stellar medium in ^''Al . Including rotation somewhat modifies 
this picture, as we will show in the following sections. 
Before entering into the details of the effects of rotation on 
^^Al release by WR stars, let us briefly recall those effects of 
rotation on massive star evolution which may be relevant in the 
context of Al nucleosynthesis : 

1 . Rotation allows diffusion of the chemical species from the 
convective core to the radiative envelope and vice versa; 
this results in changes of the surface abundances already 
during the Main Sequence phase. The surface enrichments 
given by the models are generally in good agreement with 
the observed ones (see Heger & Langer 120001 Meynet & 
Maeder l2000l Maeder & Meynet 2001] 

2. Rotation increases mass loss (see Sect 2). More precisely 
a rotating star of given initial mass and position in the HR 
diagram undergoes a higher mass loss than a non-rotating 
star (paper VI). 

3. Diffusion of elements actively contributing to the nuclear 
energy production inside the star together with enhanced 
mass loss allow massive stars to enter the WR phase at an 
earlier stage of their evolution. As a consequence, the WR 
phase lasts for longer than in the case of non-rotating stars. 
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Fig. 2. Variation of the mass fractions of -^Mg and ^''Al as 
a function of the lagrangian mass inside a rotating and non- 
rotating 60 Mq stellar model at Z =0.020 during the core H- 
burning phase. In both models the mass fraction of hydrogen in 
the convective core is equal to 0.4. 

4. Finally, rotation decreases the minimum initial mass for a 
single star to enter the WR phase (see Table 2). 

These effects, when accounted for in stellar models, much 
improve the agreement between theory and observations for 
what concerns the characteristics of the WR populations at 
diff'erent metallicities (see paper X and XI). Through these 
diff'erent processes, rotation also modifies the quantity of 
^^Al ejected by massive stars undergoing strong winds. 

In Fig. n to Fig. |5] we present the evolution of the 
^^Al abundance in test models of 60 and 120 Mq at different 
metallicities. In Fig.0]and Table 1, we provide the wind ejected 
mass of Al 
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for models with M, e [25 M© ;120 Mq ] and metallicities from 
Z - 0.004 to Z = 0.04. Here we do not take the radioactive 
decay of ^^Al into account, provided that it is the production 
rate of ^*A1 that is needed to estimate the contribution of WR 
stars to the present galactic mass of Al . 

The effects of rotation on ^^Al nucleosynthesis are the fol- 
lowings: 

1. As can be seen from Fig. ^ rotation slows down the 
decrease in mass of the convective core during central 
hydrogen burning phase. This results from the fact that 
rotational mixing continuously brings fresh H-fuel into 
the core. The regions where ^^Al is synthesized are closer 
to the surface, making it easier for the mixing to connect 



A. Palacios et al.: -*A1 production in rotating Wolf-Rayet stars 



5 



Table 1. Wind ejected mass of -^Al in units of 10 '*Mq for both rotating and non-rotating models at the different metallicities 
considered. Bold values are for models computed with metallicity dependent mass loss during the WR phase. 
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both regions, thus favouring an increase of the quantity 
of ^^Al ejected. One notes also from Fig. [l] that whether 
rotation is taken into account or not, the maximum content 
of Al inside the convective core is about the same in both 
cases, and it is reached at 1 Myr Afterwards, the amount 
of ^^Al slowly decreases as its production via the reaction 
^^Mg(p,y)^^Al is not sufficient anymore to compensate its 
reduction due to radioactive decay. 

2. Due to rotational mixing, the ^^Al surface abundance 
begins to change at a much earlier time in the rotating 
model. This can also be seen from Fig.|2lwhich shows the 
variations inside rotating and non-rotating 60 M© stellar 
models of the mass fractions of ^^Mg and ^*A1 at the 
middle of the core H-burning phase. While the abundances 
of these two elements are very similar, the main differences 
arise in the radiative envelope: in the non-rotating model, 
above the transition zone left over by the recession of 
the convective core, the abundances present a flat profile 
at a level corresponding to the initial abundance (which 
is of course for ^^Al ). Instead in the rotating model, 
the inwards diffusion of ^""Mg and the outwards diffusion 
of ^^Al smooth the profiles and produce changes of the 
surface abundances while mass loss has not yet uncovered 
layers previously pertaining to the convective core. 

3. From Fig. and Fig.|3l it appears that most of the ^^Al is 
ejected while the star is at the beginning of the WR phase. 
Indeed at this stage, the layers removed by the stellar 
winds are heavily loaded in ^^Al and the mass loss rate 
is high. When layers processed by He-burning appear at 
the surface, the quantity of ^^Al much decreases, because 
^^Al is destroyed (mainly by neutron capture) at the 
beginning of the core He-burning phase. Since rotating 
stars enter the WR phase at an earlier evolutionary stage 
than non-rotating ones (see Fig. [Q, the period during 
which ^^Al is released into the ISM is longer, and the net 
yield of this nuclide is also increased when rotation is taken 
into account. This appears clearly in Table 1, where we 
show the total wind ejected mass of ^^Al for the computed 
models. 

4. The lowering of the minimum mass for a star to enter the 
WR phase contributes to enhance considerably the total 
contribution of WR stars to the global galactic content of 



Al . Less massive stars, favoured by the IMF, will now 
contribute to the ^^Al enrichment of the medium through 
their winds, whereas their non-rotating counterparts could 
only contribute through the supernova explosion (see Sect. 
4). 

3.1. Relative effects of mass loss and rotation on ^^Al 

Values listed in Table 1 show that for the higher metallicity 
models the ^^Al mass ejected is increased by a factor of about 
2.5 for typical WR progenitors with Minn ~ 40-85 Mq . For 
the 25 Mq model, the ^^Al mass released by rotating models 
is 33 times higher than the one obtained for non-rotating ones, 
mainly because in the latter case, these stars do not enter the 
WR phase and contribute very few through their winds to the 
^^Al enrichment of the interstellar medium. In the lower mass 
stars, the wind during the WR phase is weaker than for the more 
massive ones, and rotation-induced mixing is the main mech- 
anism driving the evolution of the surface abundance patterns. 
On the contrary, in the mass range ]60 Mq ; 120 Mq ], the sur- 
face abundances during the WR phase essentially result from 
the efficiency of the wind to peel off the star, so that the effects 
of rotation-induced mixing do not appear that strikingly when 
comparing rotating and non-rotating models (see also Fig.|3i. 

It is interesting to note from Fig. |4] that the yields from 
Meynet et al. (1997), computed for non-rotating models with 
enhanced mass loss are lower than the ones obtained for rotat- 
ing models with lower mass loss rates compatible with the most 
recent observational determination (see Sect. 2). This again 
points out the importance of rotation, non solely in increasing 
mass loss, but mainly in allowing for transport of the chemicals 
throughout the stellar interiors. 



3.2. Metallicity dependence of °Al yields 

As shown in previous papers (Meynet et al. ll997l Vuissoz et al. 
i2004j . the amount of ^^Al both synthesized and expeUed by a 
WR star increases with metallicity. This is already the case for 
non-rotating models, and rotation adds to this effect as can be 
seen in Figs. 3 and 4. 

In Fig. 13 we present the evolution of the surface and cen- 
ter mass fraction of ^^Al as a function of time for rotating 120 
Mq stellar models at various metallicities. One sees that the 
higher the metallicity, the larger the amount of ^^Al produced 
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Fig. 3. Same as Fig. 1 for rotating models of a 120 Mq star at four different metallicities, as indicated on the plots. 



in the convective core and the earlier ^^Al appears at the sur- 
face. Thus for a given initial mass and rotation, the mass of 
^^Al ejected is greater at higher metallicity. This is also well 
illustrated in Fig.|3 From the present models one expects a re- 
lation of the type F(^^A1, Mi„i) oc Z^. between the yield and the 
metallicity. Here, the value of fi depends on both initial mass 
and rotation. Typically from the rotating stellar models with 
Mini = 40, 60 and 120 Mq , one obtains /J = 1.4, 1.8 and 1.4 re- 
spectively in the metallicity range from Z=0.008 to 0.040. For 
the non-rotating 60 and 120 M© and in the metallicity range 
from 0.02 to 0.04, one obtains values of fi equal to 1.2 and 1.9 
respectively. As can be seen from these numerical examples, 
the concomitant effects of rotation and mass loss are intricate 
and it is difficult to draw more general trends for what concerns 
the sensitivity of the yields on the initial metallicity. 



loss rates are multiplied by a factor ~ 1 .4. The resulting wind 
ejected masses of ^^Al are listed in bold face in Table 1 . For 
the most massive stars, the masses of ejected ^^Al do not 
differ by more than 13 % when the metallicity dependent mass 
loss rates are used. At first sight this may appear surprising. 
However the most massive stars enter the WR phase while 
central H-burning is still going on. In that case a stronger mass 
loss produces a more important reduction of the convective 
H-burning core, accompanied by a severe reduction of the 
central temperature. The end of the H-burning phase is then 
lengthened, and thus the ^^Al present inside the convective 
core has more time to decay in situ before being expelled. 
The changes are however small and therefore one can conclude 
that the use of metallicity dependent mass loss rate does not 
affect significantly the results. 



3.3. Effects of metallicity dependent mass loss during 
the WR phase 

For twice the solar metallicity, we computed models using 
Crowther et al. ti2002j prescription for the mass loss rates 
during the WR phase. For this metallicity, the derived mass 



3.4. Yields dependence on initial velocity 

Finally let us recall that in Table 1, the yields in ^*A1 for ro- 
tating models are given for an initial rotation velocity on the 
ZAMS of 300 km.s"' . This value corresponds to an average 
velocity during the MS phase between 200 and 250 km.s"' , i.e. 
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Fig. 4. ^''Al yields, as listed in Table 1, as a function of ini- 
tial stellar mass and metallicity. Filled symbols are for rotat- 
ing models with = 300 km.s open black squares and 
triangles are for non-rotating models and open grey circles 
connected by long-short dashed lines are from Meynet et al. 
(1997). The corresponding metallicity is indicated on the right 
part of the figure. 

well in the observed velocity range for the OB Main Sequence 
stars at solar metallicity. In order to estimate the effects of a 
higher initial velocity, we have computed, for solar metallicity, 
a 60 Mo stellar model with Wini - 500 km s The quantity of 
^^Al ejected by the stellar winds is in this case of 2.6 10""* M©, 
which is less than 20% higher than the yields obtained from 
the model with Uini = 300 km s Thus at least in this par- 
ticular case, the yields do not appear to be very sensitive to a 
change of the initial velocity in the range between 300 and 500 
km s"'. 

4. Contribution of WR stars to the ^''Al galactic 
content 

Let us first consider a simple way to estimate the global con- 
tribution of WR stars to the present day ^^Al content of the 
Galaxy. According to van der Hucht ( 2002 ), the galactic dis- 
tribution of WR stars shows a projected surface density of 2.7 
WR stars per kpc^ in the solar neighbourhood. Extrapolating 
this estimate of the surface density to the whole Galaxy gives 
a total number of WR stars in the Galaxy of the order of about 
2000. Since the average lifetime of WR stars is ~ 1 Myr (e.g. 
paper XI) this gives a galactic frequency /wr ~ 2000 Myr"^ 
Let us note that these numbers are likely to be lower limits, 
since the star formation rate and the number of WR to O-type 
stars is much higher in inner galactic regions. 

The galactic frequency of WR stars may also be estimated 
by using a stellar Initial Mass Function (IMF) and the current 



death rate of massive stars in the Milky Way. Observations of 
external spiral galaxies indicate that the expected supernova 



frequency in our Galaxy is /sn = 2.5_y5 per century, with 85% 
of them coming from massive stars (Tammann et al. 1994). We 
adopt then /sn - 2 per century for the death rate of massive 
stars in the Milky Way. We use a Salpeter J1955> stellar IMF-', 
<b{M) oc M with X - 1.35. The WR frequency is then 
given by: 



/^"^ <D(M)^M 

r _ r -^^LWR 

/WR - JSN 



/sn 



(2) 



where Mlsn and Mlwr are the lowest masses of stars that 
may become supernovae and WR, respectively. The former 
is Mlsn ~ 8 M©, while the latter depends on metallicity and 
adopted mass loss rates. According to our calculations with ro- 
tation Mlwr ~ 21 Mq for Z = 2Zq (see Table 2). We consider 
here the high metallicity models since, most of the WR stars 
are found in the inner regions of our Galaxy where the metal- 
Ucity is higher (see below). Taking these values into account 
one obtains /wr ~ 0.2 /sn, i. e. a WR frequency of 0.4 per 
century or 4000 WR stars per Myr. On the other hand, from the 
yields of Table 1 and our adopted IMF, an average WR yield of 
<M^) ~ 3.4 lO^'^Mo can be defined for Z = 2Zo stars. 

The collective production rate of Al in the Galaxy by WR 
stars can then be expressed as 



M26 ~ 1.4 



/v. 



3.410-4 4000 



Mq ■ Myr' 



(3) 



In the case of a stationary regime as it assumed to be the case 
in the Milky Way, the galactic mass of Al does not vary with 
time, so that at a given time t, the production {Pieif)) and de- 
struction (D26(f)) rates of ^^Al are equal. The destruction rate 
is easy to evaluate, considering the evolution of the ^^Al mass 
solely due to its decay , e.g. 



M26(0 = M26(fo) ■ e- 



£>26 = 



M26 



As ^26(0 = D2(,(t), the production rate of Al is then given by 



^26 ~ ^26 ~ '"26 



\l 



gal 



Mq per unit time. 



(4) 



The total mass of Al originating from WR stars present in the 
Galaxy is related to the production rate of ^^Al by these stars 
through Eq. (4). The mean Ufetime of ^^Al being of the order 
of 1 .04 Myr, we finally have: 



M?f = T26XMf 



^26 



'26 



(5) 



Thus Eq. (3) (and Eq. (6) in the following) in addition to give 
the injection rate of ^*A1 by WR stars, can also be used to 
evaluate the total mass of ^*'A1 originating from WR stars. 
According to our rough estimate from Eq. (3), the mass of 



^ In the following, we will also discuss the effect of using a larger 
slope (x = 1.7) that could be more appropriate for massive stars when 
accounting for binarity corrections (Scalo 1986; Kroupa et al. 1993). 
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Table 2. Minimum mass Mmrj^ for a star to enter the WR phase 
as a function of Z and y,„,- (from Meynet &■ Maeder (2004)) 
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Z 


Z 




(km.s-' ) 


0.04 


0.02 


0.008 


0.004 


Mlwr 





29 


37 


42 


52 


(Mo) 


300 


21 


22 


25 


32 



^^Al contributed by WR stars amounts to 45 % and even 
70 % of the observationally inferred mass of about 1.8 
(Knodlseder il997ll l. 

In a more formal way, the ^^Al production rate (in solar 
mass per unit time) from galactic WR stars may be evaluated 
as 

M26 = I 2nro-(r) I (^(M)MY(,[Mi,Zi(r)]dM dr, (6) 
Jo Jmlwr('-) 

where cr(r) is the Star Formation Rate (hereafter SFR) per unit 
surface in the Galaxy given as a function of galactocentric dis- 
tance r and is the mass of ^^Al ejected by a star of initial 
mass Mi and initial metallicity Z,(r), according to Eq. (1). 
To account for the metallicity dependence of the ^^Al yields 
and of the lowest WR mass Miwr in Eq. (6), we use the Z- 
dependent values given in Table 1 and in Table 2 respectively. 

The integration in Eq. (6) requires the following ingredi- 
ents: 

- • Normalization of both the SFR cr(r) and the IMF <I)(M). 
We can either rely, as in our first simple estimate, on an ap- 
proximate frequency of SNae events in the Galaxy; in this 
case, the integral 2nrcr(r)dr <i>(M)dM is normal- 
ized to 2 10"^ Myr ' (2 cc-SNae per century in the Galaxy). 
We can also follow Knodlseder ( 1999 1, and normalize the 
SFR so that the total predicted Lyman continuum pho- 
ton production of the Galactic population defined by the 
IMF matches the observed Lyman continuum luminosity 
Q. The actual value of Q as deduced by McKee & Williams 
<1997> is equal to 2.56 xlO^^ photons, s"'. In that case. 



L 2nro-(r)dr normalized to --po- 



according to Eqs. (18) and (20) of Knodlseder J1999> . we 
will have (b(M)dM - 1 and the star formation rate 

— , sh(M) 

/g " SH{MmM)T,„,AM)dM 

being the time-averaged Lyman continuum luminosity and 
TtotaKM) the total lifetime of a star of initial mass M. To 
evaluate these two quantities, we make use of Eqs. (21) 
of Knodlseder ( 11999]) . which is a rough approximation 
when taking rotation into account, considering that these 
estimates are for non-rotating models with diff'erent mass 
loss prescriptions than the ones used in the present paper 
Keeping this in mind, we nevertheless found it interesting 
to estimate the effect of normalization on the ^^Al injection 
rates predicted by Eq. (6), and used both SNae and Lyman 
continuum for this purpose. 
- • The metallicity gradient plays an important role for de- 
termining the global contribution of WR stars to the galac- 
tic ^^Al enrichment. In order to illustrate quantitatively this 



point, we performed three estimates with three different 
galactic metallicity distributions in the Galaxy. 
Our first choice (model Zl) is based on the metallicity gra- 
dient deduced from observations of HII regions by Shaver 
et al. (I1983> and on the galactic centre metal abundance 
obtained by Najarro (I1999i . It corresponds to a metallic- 
ity gradient of d log Z/dr - -0.07 dex.kpc"' (with Z(8.5 
kpc)=0.02) in the whole region lying between 1.7 and 15 
kpc. In the inner 1.7 kpc, the metallicity was taken equal to 
0.06 (i.e. three times solar). Since we did not compute mod- 
els at such a high metallicity, we estimated the ^^Al yields 
for metallicities above 0.04 extrapolating those obtained at 
solar and twice the solar metallicity. 

Our second choice (model Z2) is based on the metallicity 
gradient proposed by Andrievsky et al. ( 2004 1, which relies 
on cepheids observation. They divide the Galaxy into three 
zones: Zone I encompasses the region between 4.0 and 6.6 
kpc. In this zone metallicity varies according to the law 
log Z/Zq = -0.128Rc + 0.930, Rg being the galactocentric 
distance expressed in kpc. Zone II extends from 6.6 to 10.6 
kpc and has a metallicity law log Z/Zq - -0.0447;c + 0.363 
(for Re - 8.5 kpc, Z ~ Zq). The third zone goes from 
10.6 to 14.6 kpc and has a metallicity law of the form^ 
log(Z + 0.024)/Zq = +0.0047;c + 0.256. For galactocen- 
tric distances inferior to 4 kpc, we take a constant value of 
the metallicity equal to 2.62 x Zq. 

Finally in model Z3, we consider a much shallower gra- 
dient of dlogZ/dr - -0.04 dex.kpc ' in the whole re- 
gion between 4.4 and 12.9 kpc as suggested by Daflon & 
Cunha (2003 1. For galactocentric distances below 4.4, we 
suppose a constant metallicity equal to 0.03 and for galac- 
tocentric distance above 12.9 kpc, we take a constant metal- 
licity equal to 0.013. 

These three hypotheses not only account for the various 
values of the metallicity gradients found in the literature, 
but also for the different views on the metallicity of the 
galactic centre. For instance Ramirez et al. ( I2000> and 
Najarro et al. (p004i find that the mean [Fe/H] of the 
Galactic Centre stars is similar to that of the solar neigh- 
bourhood stars, while Najarro ( II999> obtains values be- 
tween 2 and 3 times the solar metallicity. 
• Different estimates can be found in the literature for the 
radial distribution of the present Star Formation Rate in the 
Galaxy cr(r). We chose to consider two different prescrip- 
tions : Paladini et al. ti2004> . give the distribution of HII 
regions in our Galaxy as a function of the galactocentric 
distance. HII regions are considered to be good tracers of 
the SFR, and we suppose that their surface density distribu- 
tion is directly proportional to it. Wang & Silk ( 1 994 1, on an 
other hand, give the SFR per unit surface, based on obser- 
vations of various tracers (Lyman continuum photons from 
HII regions; pulsars; supernova remnants) and normalized 
to the solar neighbourhood. 



There is an error in the expression given by Andrievsky et al. 
<2004> for this zone, that can be corrected considering log(Z + 
0.024) /Zo instead of log Z/Zq. 
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Applying Eq. (6) with a Salpeter IMF slope, a SFR from 
Wang & Silk ( I1994> and a normalization on the rate of SNae 
events, leads to an injection rate of M ~ 0.3 - 0.6 Mq /Myr of 
Al for the non-rotating models and to values between 0.9 and 
1.3 Mq /Myr for the rotating models, not far from the order of 
magnitude estimates performed above (see Eq. 3). Rotation, all 
other physical ingredients being the same, thus roughly dou- 
bles the contribution of the WR stellar winds to the galactic 
^^Al enrichment, as can be seen from values in Table 3. This 
remark holds whatever the combination of the physical ingre- 
dients entering Eq. (6) used. In that respect, rotation (for the 
adopted value for the initial velocity, Uini - 300 km.s"' ) has a 
similar impact on the global contribution of the WR stars as an 
enhancement of mass loss rates (see Meynet et al. ll997l l. 

In Fig.|5j Fig.|6l Fig. and Fig.|S] we compare the theo- 
retical mass surface density profiles of ^^Al , E26 (dashed his- 
tograms) with the profile deduced from COMPTEL data anal- 
ysis (Knodlseder l 1 9971 black curve). In Fig.|5] we use the SFR 
from Wang & Silk (1994), normalized to a rate of 2 SNae per 
century. In Fig.|6land Fig.0, the SFR is estimated following 
Fig. 3 of Paladini et al. J2004L and it is normalized either on 
the SNae rate (Fig. |6|l, or on the total galactic Lyman contin- 
uum luminosity (Fig.Qi. For each case we compare the effect 
of rotation and of the adopted galactic metallicity gradient. In 
Figures 5 to 7, we used a Salpeter IMF; Fig.|S]is analogous to 
Fig -El (right panels) and Fig. (right panels) but for a Kroupa 
IMF ((D(M) oc M-2-^). 

- The theoretical predictions are in general below the ob- 
served surface density curve. This is expected since in our 
estimate we only account for the contribution of the winds 
of WR stars, while other sources as the supernovae likely 
also contribute. 

- In the present theoretical framework, the fraction con- 
tributed by other sources could be estimated from the dis- 
crepancy between the theoretical and observed histograms. 
From Fig. |5] to Fig. we see that the contribution of 
other sources than WR stars is much larger when the yields 
of non-rotating models are accounted for. As recalled in 
Sect. 1 (see paper XI for more details), the non-rotating 
models fail to fit many observed properties of WR star pop- 
ulations that rotating models succeed in reproducing. We 
will thus preferentially trust the results obtained from the 
rotating models. 

- The higher the metallicity of the inner galactic regions 
(r < 6), the more important the WR contribution to the 
^^Al synthesis. In particular, models labelled as Z3, with 
the shallow metallicity gradient and an adopted value of 
0.03 for the metallicity in the innermost regions, lead to the 
lower ^^Al injection rates in all the cases presented here. 

- In the outer regions, beyond 8 kpc, there is a deficit in all the 
models considered here when the initial rotation velocity 
assumed is that of Viai - 300 km s"^ . If the initial distribu- 
tion of rotation favours faster rotation at lower metallicity, 
such a deficit could be smaller, but at the present time the 
observational evidence for such a behaviour remains indi- 
rect (Maeder et al. .l999) . 



Salpeter : $(M) « 



SFR Wang k Silk 
Model Zl 
v.^. = 300 km s- 
l!f= 1.242 M, 




2 4 6 8 10 2 4 6 8 10 12 
Galactocenlric radius (kpc) Galactocentric radius (kpc) 



Fig. 5. Mass surface density profile of Al (in 10" ■ M© -kpc" 
) as a function of the galactocentric radius for an IMF index 
of 1.35 and a SFR from Wang & Silk (1994), normahzed in 
order to reproduce a rate of 2 SNae events per century in the 
Galaxy. The black histogram with error bars is the profile de- 
duced from COMPTEL data analysis (Knodlseder 1997). The 
dashed histograms are the predicted distribution for the pre- 
scriptions adopted when using non-rotating (left panels) or ro- 
tating models (right panels). 



- Around 50% of the total mass of ^^Al in the Galaxy origi- 
nates from the regions between 3 kpc and 6 kpc, whatever 
the metallicity gradient adopted. This is related to the shape 
of the SFR, which presents a peak in this region (see Fig. 3 
in Prantzos & Aubert 119951 and in Paladini et al. I2004> . 
mainly due to the presence of a ring of molecular clouds 
around 5 kpc. 

- The important dependence of the theoretical results on 
the adopted metallicity distribution in the Galaxy appears 
clearly through these numerical experiments. Depending on 
the metallicity distribution law adopted, both the radial dis- 
tribution of the Al and the total quantity of ^^Al predicted 
to be deposited by the WR stellar winds vary significantly. 
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Fig. 6. Same as Fig. 5, for a SFR based on Paladini et al. (2004) 

rR rl20 

number counts of HII clouds, and 2nro-(r)dr (S>(M)dM 
normalized to 2 lO'* SNae/Myr. 

Shallower metallicity gradients, as adopted in model Z3, 
produce smaller values for and smoother profiles of 
^26. This conclusion holds whether rotation is included or 
not. 

Thus from the observed distribution of S26, it is not pos- 
sible to infer constraints on the metallicity dependence of 
the stellar yields (even if only one source dominates), un- 
less more reliable and concordant estimates of the galactic 
metallicity gradients are obtained. 
- Comparing Fig. |5] Fig. |6] and Fig. also allows us to 
study the effects of the choice of the SFR or that of its 
normalization on the ^^Al distribution and injection rate by 
WR stars. The SFR adopted entirely determines the shape 
of the surface density distribution predicted : while the 
Wang & Silk (1994) prescription leads to a monotonous 
decrease of E26 with increasing galactocentric radius, the 
values derived from Paladini et al. ( 2004 1 all present a low 
surface density in the inner regions, a peak in the region 



Fig. 7. Same as Fig. 5, for a SFR based on Paladini et al. 
(2004) number counts of Hll clouds, normalized to the ob- 
served Lyman continuum luminosity of Q = 2.56 lO^"* photons 

of the molecular clouds ring, and then a monotonous 
decrease, in better agreement with the nominal observed 
distribution^. On the other hand, as can be seen from Table 
3, these differences between the two prescriptions used for 
the SFR (Fig. |5l and Fig. |6j only slightly affect the total 
contributed mass of ^''Al predicted, both for rotating and 
non-rotating WR models. 

Changing the normalization, and adjusting the SFR to the 
match the observed galactic Lyman continuum luminosity 
does not change the results significantly, despite the fact 
that the value adopted of 2.56 photons .s ' coiTesponds 
to a lower rate of SNae events, namely 1.2 SNae per 
century. The masses deduced with the same value for 
Q as Knodlseder (1999), i.e. 3.5 photons .s"' are to be 

' Note that within the errorbars from observations, the shape ob- 
tained with the Wang & Silk lil994il prescription also matches the ob- 
served distribution. 



A. Palacios et al.: -*A1 production in rotating Wolf-Rayet stars 



11 



Kroupa : ^(M) « M "^^ 



1.5 



"0.5 




"0.5 



-JTTTy 



SNae normalization 
lifodel Zl 

t),j, = 300 l<m S-' 
0,759 M, 



SNae normalization 
Model Z2 
ti;„ = 300 l<m S-' 
te'= 0,838 M, 



SNae normalization 
Model Z3 

= 300 l(m s-i 
M!J'= 0,565 M„ 




Lyc normalization 
Model Z2 
ti^„i = 300 km S-' 
M!j= 1.195 M„ 




Lyc normalization 
Model Z3 
D,„i = 300 km s-i 
Mg= 0.805 Mg 



2 4 6 8 10 2 4 6 8 10 12 
Galactocentric radius (kpc) Galactocentric radius (kpc) 



Fig. 8. Same as right column of Fig. 6 and Fig. 7 for a Kroupa 
IMF 



multiplied by 1.37, which leads in the most favorable case 
to 1.70 M0 of the present galactic mass of ^^Al originating 
from WR stars. 

Finally, using an IMF slope of x = 1.7 (Scalo, 1986; 
Kroupa 1993) rather than 1.35 lowers the total galactic 
mass of ^^Al predicted by the models as was shown in 
Fig.|S]and in Palacios et al. ( (20(J 4lin their Fig. 3), without 
affecting the shape of the ^^Al surface density distribution. 



From the above estimates and comparisons, we can con- 
clude that when stellar models accounting for the observed 
properties of the WR star populations are used (i.e. rotating 
models), WR stars appear to be significant if not dominant con- 
tributors to the present day ^^Al content of the Milky Way. The 
predictions for the total mass of ^^Al originating from WR stars 
vary from 0.6 Mq to 1 .4 Mq (see Table 3), whereas for the black 
histogram to which we compare our predictions, Knodlseder 
(1997) obtained a total galactic mass of ^''Al of 1.8 + 0.1 Mq . 
The values listed in Table 3 show that WR stars can account for 



Table 3. ^^Al injection rates predicted by the different mod- 
els computed : (A) Salpeter's IMF, SFR from Wang & Silk 
(1994), SNae normahzation (Fig. 5); (B) Salpeter's IMF, SFR 
from Paladini et al. (2004), SNae normalization (Fig. 6); (C) 
Salpeter's IMF, SFR from Paladini et al. (2004), Lyc normal- 
ization (Fig. 7); (D) Kroupa's IMF, SFR from Paladini et al. 
(2004), Lyc normalization (Fig. 8); (E) Kroupa's IMF SFR 
from Paladini et al. (2004), SNae normalization (Fig. 8). 
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20% to 93% of the observational estimates of the actual total 
galactic ^*A1 mass (1.5-3 Mq ). The predicted masses are little 
affected (w 10% variation) by the uncertainties on the galactic 
SFR, the IMF slope or the quantity against which the number 
of massive stars is to be normalized. On the other hand, the pre- 
dicted total mass of ^^Al originating from WR stars varies by 
as much as 30% according to the galactic metalhcity gradient 
adopted. The shape of the predicted surface density distribu- 
tions of ^^Al , S26, against the galactocentric radius strongly 
depends on the SFR adopted, so that the comparison with the 
values deduced from observations is of little interest. 



5. Conclusion 

We have studied the impact of rotation and initial metalhc- 
ity on the yields of ^^Al ejected by the WR stellar winds. 
When making use of rotation velocities corresponding to the 
observed average rotational velocities during the MS phase, 
i.e. for values around 200 km.s ' , the inclusion of rotation, 
all other things being equal, roughly doubles the WR contribu- 
tion. Higher initial rotation velocities would still increase the 
quantity of ^^Al ejected by the WR stellar winds. 

We compared the effects of rotation and those induced by 
an enhancement of the mass loss rates. On one hand, non- 
rotating models with enhanced mass loss rates (Meynet et al. 
il997ii as well as the rotating models presented here (for which 
the mass loss rates are quite small), are both able to reproduce 
the observed variation with the metallicity of the number of 
WR to O-type stars. On the other hand, the ^^Al yields pre- 
dicted are different, especially in the lower mass domain. This 
is evidence of the fact that rotation has other more subtle effects 
rather than an enhancement of the mass loss rates, as already 
pointed out in Sect. 3. Our rotating models thus predict wind 
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ejected masses that are about a factor of 1.4 superior to the 
ones provided by Meynet et al. (1997). 

We showed that the use of metallicity dependent mass loss 
rates during the WR phase as suggested by Crowther et al. 
(2000) does not have a significant impact on the results. 

The global contribution of WR stars remains difficult to as- 
sess in view of the numerous uncertainties pertaining not only 
to some physical ingredients of the models as the mass loss 
rates or the nuclear reaction rates, but also to galactic param- 
eters such as the galactic metallicity gradient or the supernova 
rate. We showed that the choice of the SFR tracer determines 
the shape of the ^^Al surface density distribution with galac- 
tocentric distance, but marginally affects the predicted con- 
tributed mass of ^^Al by WR. On the other hand, the choice 
of the metallicity gradient as well as that of the IMF slope 
may significantly aflFect the contribution of WR stars to the total 
galactic mass of ^*A1 , which can amount to 0.6 (shallow 
metallicity gradient with Kroupa's IMF) up to 1.4 M© (steep 
metallicity gradient with Salpeter's IMF). 

In that respect the study of associations such as the Cygnus 
region, sufficiently young to prevent supernovae having much 
contributed to the ^^Al enrichment are of prime importance. 
Possibly the detection of the 1.809 MeV emission around a 
single source would be even more important. In that respect 
the upper limit found for y^-Vel (Oberlack et al. 2000) already 
provides a strong constraint. The binary nature of the object 
however requires some care when comparisons are made with 
a single star model. This particular case will be studied in detail 
in a forthcoming paper. 
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